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BACKGROUND OF THE INVENTION 

Field of the Invention 

[0002] The present invention relates to apparatus for monitoring a system with time in 
space and method therefor and more particularly relates to novel apparatus such as a 
unique MRI machine, a novel irrigation apparatus for testing the effectiveness of an 
irrigation system, a novel testing apparatus for determining the efficiency of a heating 
and/or cooling system, and the like, for testing or controlling a system in which fluid 
flows and where in the system fluid either dissipates in part or requires regeneration. 

Description of the Prior Art 

[0003] Presently apparatuses are known for monitoring testing or measuring a system 
in which a fluid that is flowing or substances in the fluid will dissipate in part as it 
traverses the system or will require regeneration. For example. MRI machines are 
used today to create images with or without administration of a tracer-contrast agent. 
Customarily, the machine is controlled to take a series of images at discrete time 
intervals and the images are then dynamically analyzed to obtain an output result. For 
example, dynamic studies of contrast enhancement in breast tumors have 
demonstrated that the rate of change in signal intensity is an important parameter for 
the distinction of breast masses, leading to pharmacokinetic studies. However, it is 
known that as a result of tumor heterogeneity, there are significant local variations in 
the time evolution of contrast enhancement, and, therefore, maintaining high spatial 
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resolution in both the recording and analysis steps is very important. In a standard 
clinical MRI of the breast, it is difficult to achieve high spatial resolution and also 
maintain high temporal resolution. In most dynamic studies performed previously, the 
emphasis was on high temporal resolution (at the expense of spatial resolution) 
monitoring the equilibration in the intravascular space and early diffusion into the 
extracellular space of the tissue. As a consequence, in standard MRI machines the 
output results are sometimes inconclusive. The foregoing is also characteristic of other 
systems in which a fluid flows or a component thereof dissipates in part or requires 
regeneration, such as, for example, an irrigation system, a heating and cooling system 
and the like. 

SUMMARY OF THE INVENTION 

[0004] Accordingly, the object of the present invention is to provide an apparatus or a 
machine, and a correlated method, for monitoring a system, in which a fluid is flowing, 
with time in space, which will provide more conclusive results regarding system 
anomalies or system efficiency. 

[0005] The present invention relates to an apparatus for monitoring a system with time 
in space. The system can be physical, chemical, biological, physiological, 
environmental, clinical or any other system in part or in whole, the system evolving with 
time over space in a certain way. The apparatus of the present invention can function 
on the basis of one, two, three or higher dimensions. The type and extent of spatial 
resolution and the number of time points and their spacing, that the apparatus selects, 
depend on the system and can be varied with a lower limit for the number of time 
points of two. For example, it can be used for processing time dependent data of 
radiologic examinations such as MRI, ultra-sonography. X-ray tomography or 
conventional X-ray, or Nuclear medicine for obtaining diagnosis, prognosis and therapy 
follow up of tumors or any other pathological disorders. It can be utilized for processing 
monitoring or controlling environmental data of water irrigation. It can be used to 
analyze data that will permit determination of leafing areas in pipes. It can be used to 
analyze data obtained in the food, cosmetic and other industries which involve mixture 
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and solution preparations and detemiination of their homogeneity. It can be also used 
to assess the efficiency of heating and/or cooling systems. 

[0006] There are numerous phenomena that evolve over space with time in a way that 
can be treated according to the present invention by utilizing a novel approach which is 
termed herein as by wash-in and wash-out behavior The wash-in and wash-out are 
terms that are used symbolically to describe a change in one direction (wash-in) and 
the reverse change (wash-out) which may not be true reversal but can be any pathway 
that induces a change. Specifically flow of fluid in a system where the fluid or fluid 
component dissipates or needs to be regenerated, is described according to the 
invention as wash-in and washout. 

[0007] For any wash-in/wash-out situation, it is possible to describe n numbers of 
patterns of wash-out, when n can range from 1 to any integral number, 2, 3, 4 etc., on 
the basis of m time points, when m can range from 2 to any integral number of specific 
time points in the time evolution of the process. The definition of wash-out is not strict 
and a wash-in can become a wash-out and vice versa. 

[0008] For each kind of system, the apparatus of the present invention provides means 
for monitoring, controlling or regulating the system by providing means for setting time 
points and other optimal parameters of the system. This setting uses a novel 
calibration map based on a physical model which describes the evolution with time in 
an approximate or rigorous manner. These calibration maps serve also to interpret 
quantitatively the final color hue/color intensity coded maps obtained as one of the 
products of the apparatus. 

[0009] One particular use of the novel apparatus is for contrast enhanced MRI data in 
order to obtain products that facilitate specific diagnosis of cancer. The time of start of 
contrast administration is time point tO and then two post contrast time points tl and t2 
are utilized. These post contrast times are selected by constructing calibration maps 
based on modeling the kinetics of contrast enhancement that relates the 
wash-in/wash-out rates to two pathophysiological parameters: microvascular 
permeability times surface area (termed in short, microvascular pemneability and 
represented by the letter K) and fraction of extracellular volume represented by v. The 
calibration map is constructed by the apparatus. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0010] Fig. 1 is a block diagram of tlie novel apparatus of the present invention for 
monitoring a system with time in space; 

[0011] Fig. 2 is a block diagram of the details of the selecting means of the apparatus 
in Fig. 1 for setting the parameters for data collection; 

[0012] Fig. 3 is a block diagram of the details showing the colorize function C of the 
selecting means of the apparatus of Fig. 1; 

[0013] Figs. 4a and 4b are a flow diagram of the selecting means of the apparatus of 
Fig. 1 illustrating the setting of the parameters for data collection; 

[0014] Figs. Sa and 5b are a flow diagram of the selecting means showing the novel 
apparatus selecting means used as a novel MR! apparatus for collection of MR! 
images. 

[0015] Fig. 6 is a block diagram of the apparatus of Fig. 1 showing in detail the control 
and effect means of the apparatus for controlling and effecting data collection on a 
system; 

[0016] Fig. 7 is a block diagram of the apparatus as applied to MRI and shows the 
details of the control and effect means of the apparatus for controlling and effecting 
MRI data collection on a subject; 

[0017] Fig. 8 is a block diagram of the apparatus of Fig. 1 showing details of the 
processing means of the apparatus for processing the collected data; 

[0018] Fig. 9 is a flow diagram of the apparatus of Fig. 1 showing details of the 
processing means of the apparatus for processing data collected in two dimensions; 
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[0019] Fig. 10 is a flow diagram of the apparatus of Fig. 1 showing the details of the 
processing means used for MR! for processing images collected from a subject; 

[0020] Fig. 1 1 is a block diagram of the apparatus of Fig. 1 showing the details of the 
analysis means of the apparatus for analyzing the processed data; 

[0021] Fig. 12 is a block diagram of the apparatus of Fig. 1 showing the details of the 
analysis means used for MRI for analyzing 3TP images; 

[0022] Fig. 13 is a graph showing the pattern of slow wash-out/slow wash-in and slow 
wash-out/fast wash-in for the three time points tO, t1 and t2 and with data values I(t0), 
I(t1),l{t2); 

[0023] Fig. 14 is a graph like Fig. 13 showing the pattern of moderate wash-out/slow 
wash-in and of moderate wash-out/fast-wash-in; 

[0024] Fig. 15 is a graph like Fig. 13 showing the pattern of fast wash-out/slow washin 
and of fast wash-out/fast-wash-in; 

[0025] Fig. 16 is a schematic illustration of a typical calibration map as used in the 
novel apparatus for MRI. 

[0026] Figures 17 (a), (b) and (c) show, respectively, 

(a) MLO mammographic projection of a 49 yo woman does not show 15 
mm cancer which was palpable and diagnosed by in-office FNA as 
infiltrating ductal cancer. 

(b) Subtraction sagital MRI (0 minute (pre-contrast) image subtracted from 
2 minute image shows ring-enhancing, spiculated malignancy in the 
anterior breast (arrow). Posteriorly, a small intramammary lymph note 
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also enhances, but has a characteristic morphology, including fatty 
hilus (arrowhead). 

(c) 3TP parametric similar to that of 1 b, but overlay of colored pixels is 
superimposed on MRI image by 3TP software based on a 
physiological model (described in text). Bright red indicates high 
probability for malignancy and this lesion was prospectively given a 
score of 5 (highly suspicious for malignancy). 

(d) Calibration map corresponding to patient shown in 1a. Color hue (red, 
green, blue) is based on differences in signal intensity between the 
second and third images of the three image 3TP image set. Color 
intensity is based on the difference in SI between the first and second 
images of the image set. As described in the text, areas of high vessel 
permeability x surface area and low extravascular volume fraction 
(EVP), typical of malignancy, will be coded as bright red. For optimal 
discrimination of benign and malignant lesions, the imaging 
parameters are chosen to approximately divide the calibration map into 
equal areas of red, green and blue. 

[0027] Figures 18 (a) and (b) show, respectively, 

(a) MLO mammographic projection showing an ovoid focal lesion along 
the inferior mid breast (arrow) in a 69 yo woman. Pathological 
diagnosis, via excisional biopsy, was benign breast tissue. 

(b) 3TP sagital MRI parametric image shows virtually all dark blue pixels 
indicating a lesion with low values of vessel permeability x surface area 
and extravascular space - indicators of benignity. The prospective 
3TP diagnosis was benign (suspicion level 2) 

[0028] Figures 19 (a), (b) and (c) show, respectively, 

(a) MLO mammographic projection showing an area of ill-defined density 
and architectural distortion (arrow) in a 64 yo woman. Diagnosis, via 
excisional biopsy, was infiltrating ductal cancer. 
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(b) Subtraction MRI (image at 0 minutes (pre-Gd contrast) subtracted from 
image at time 6 minutes) showing an irregular area of enhancement 
corresponding to the mammographic lesion. 

(c) Top three images are magnified views from adjacent slices of the 
central portion of the 3TP parametric image, calculated using 0, 2 and 
6 minute MR images. The results show a visually indeterminate 
number of red pixels. The bottom three images are these same image 
locations, but 3TP images were recalculated using 0, 4 and 8 minute 
MR images. There is a shift toward an increasing number of red 
pixels, indicating malignancy to be more probable than benignancy. 
The prospective 3TP diagnosis was malignant, at suspicion level 4. 

[0029] Figure 20 shows a ROC curve derived from 3TP data in which lesion location 
was supplied to the researcher, but no other clinical information was supplied. The 
researcher was then asked to supply a number from 1 to 5 indicating probability for 
malignancy or benignity. 

[0030] Figures 21 (a), (b) and (c) show, respectively, 

(a) Optical close-up of mammographic magnification CO view - 66 yo 
woman with cluster of microcalcifications in the superior mid breast. 
Pathological diagnosis was DCIS, low to intermediate grade. 

(b) Optical close-up of subtraction sagital MRI image - 6 minutes after 
injection of Gadodiamide. Small irregular focus of enhancement 
(arrow) corresponds to cluster of microcalcifications. 

(c) Optical close-up of sagital 3TP image shows an area of predominantly 
bright green pixels thought to be benign (score = 2). A total of three 
out of eight DCIS lesions were misdiagnosed as benign by the 3TP 
method. An area of future investigation is to determine whether there 
is a specific pattern for DCIS that will allow increased accuracy of 
diagnosis. 
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[0031] Figures 22 (a), (b) and (c) show, respectively, 

(a) Ultrasound examination of a mammographically occult palpable lesion 
in a 45 yo wonnan showing a gently lobulated mass, without acoustic 
shadowing, typically of fibroadenoma, which was confinned by 
ultrasound-guided FNA. 

(b) Subtraction MIP MRI, 6 minutes after Gadodiamide injection, shows 
that the fibroadenoma confirmed by ultrasound is the largest of multiple 
enhancing smaller lesions. At the workstation, many of these smaller 
lesions showed an enhancement profile similar to that of the larger, 
palpable and sonographically confirmed lesion. In clinical practice 
evaluation of the multiple other enhancing lesions by manual 
placement of an ROI is Impractical. Though internal septations are 
said to be an important MRI sign of fibroadenoma, they were not noted 
in this patient. 

(c) 3TP parametric image shows predominantly central bright green and 
peripheral blue pixels, consistent with a benign lesion (score = 2). The 
3TP parametric image showed no other suspicious lesions, however 
due to the multiplicity of similar lesions, the patient is simply being 
followed. Confirmation of multiple benign lesions will be assumed if no 
malignancy is diagnosed after two years of mammographic follow-up. 

[0032] Figures 23 (a), (b) and (c) show, respectively, 

(a) Optical close-up of mammography (MLO projection) in a 44 yo woman 
showed a spiculated mass (straight arrow). An unexpected, 
mammographically occult, 2nd lesion was detected at the location 
indicated by the curved arrow, as a result of this patient volunteering 
for the 3TP clinical trial. Pathological diagnosis in both cases was 
infiltrating ductal cancer. 

(b) 3TP parametric sagital plane image of lesion suspected to be 
malignant by mammography. High predominance of bright red pixels 
indicates high value of the product of vascular permeability x surface 
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area, and low extravascular volume fraction, indicating high probability 
of malignancy (score = 5). 
(c) 3TP parametric image of a second adjacent sagital slice showing a 
second site very suspicious for malignancy (score 5). The radiologist 
discussed the scan with the surgeon and both sites were biopsied at 
the time of surgery confirming unsuspected multifocal malignancy. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

[0033] Referring now to Figure 1 , there is shown a block diagram of an apparatus for 
monitoring a system with time in space which embodies the present invention. This 
apparatus includes selecting means 1 for setting the time points for data collection, and 
for setting the parameters of the portions of the apparatus used for data collection, and 
for processing a colorized calibration map for analysis of processed data that relates 
variables of the system to unique wash-out/wash-in behaviors which use color hue 
.and color intensity for coding these behaviors. The selecting means is linked to a 
storage means 2 for storage the selected times, parameters and calibration map. The 
storage means 2 is linked to a control and effect means 3 for controlling the time and 
parameters of data collection determined by the selecting means 1 and for effecting 
the data collection. The means for effecting the data collection 3 are known in the art 
for which the apparatus will be used. 

[0034] The control and effect means 3 is linked to a storage means 4 for storing the 
data collected with time in space. The storage means 4 is linked to a processing 
means 5 for processing the stored collected data according to a 

[0035] novel unique wash-out/wash-in scheme which uses color hue and color intensity 
and which is programmed into the processing means 5. The processing means 5 is 
also linked to storage 2 to receive as part of its input parameters and values set in 1 . 
The processing means 5 is linked to a storage means 6 for storing the processed data. 
The storage 6 is linked to analysis means 7 for analyzing the processed data. This 
means 7 Is also linked to storage 2 for analyzing the processed data on the basis of 
the stored calibration map. 
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[0036] All the means can operate sequentially using all storage means, part of the 
storage means or none. Instead of a storage means, a direct output to input link 
between blocks 1 to 3 or blocks 3 to 5 or blocks 5 to 7 can exist. The portions of the 
apparatus need not function all at the same time nor at the same location. Referring 
now to Figure 2, there is shown a more detailed block diagram of the selected means, 
shown in block 1 of Fig. 1 . The function of the selecting means is to select time points 
for data collection, parameters of data collection and for producing a novel calibration 
map that relates washout/wash-in characteristics coded in color hue/color intensity to 
variables of the specific system being monitored. 

[0037] The inputs of the program are blocks 10, 1 1 , 12 and 13 and include the 
minimum and maximum values for variables of the system being monitored (block 10), 
input of constant parameters used for data collection (block 1 1 ), parameters used for 
data collection that can be varied (block 12) and three time points selected initially as 
an intuitive guess from general knowledge of the system as the starting point of an 
iterative procedure (block 13). The input of block 10 enables the apparatus to divide 
the system space defined by the variables into a grid, in block 14. For 1 variable, each 
grid point is a unit length. For 2 variables, each grid point defines a pixel, for 3 
variables each grid point defines a volume (voxel). The grid developed in block 14 and 
the inputs of blocks 11,12 and 13 serve to determine in block 15, for each grid point, a 
value for the parameter to be measured at preselected times termed 1(tO), I(t1 ) and 
I(t2), respectively. This determination or calculation uses an equation specific to the 
particular system being monitored that can estimate exactly or approximately the 
change with time in the value of the parameter to be measured for each set of 
variables of the system. Such equations are well known to those skilled in the art for 
particular systems. The values of I(t0) and I(t1 ) calculated in block 15 are used to 
determine or calculate in block 16, for each grid point, an intensity function according 

to the equation -^(^^) -^(^Q) j^jg intensity function represents a wash-in 

^ [I(tO)]itO-tl) 

characteristic of initial-rate of wash-in. Then in block 17, by a conventional loop, the 
grid point with a maximum intensity function is determined and is then outputted to 
block 18 as maximum intensity. In block 18, all other intensity functions in the 
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remaining grid points are normalized in reference to the maximum intensity. Next, the 
apparatus in block 19 colorizes each grid point according to the program shown in Fig. 
3. Referring to Fig. 3, function C, which determines the pattern of wash-out/color hue, 
is calculated or determined in blocks 22 and 23 in terms of I(t1 ) and I(t2), for each grid 
point, and a wash-out pattern/color 

[0038] hue is chosen according to: C = ^^^^^"^^^^^ . 

I{t\) 

[0039] Figures 13, 14 and 15 describe systematically the principles of choosing a 
wash-out pattern which is coded in color hue and of choosing wash in initial rate 
(apparent initial rate) coded in color intensity. The first pattern (Fig. 13) is the slow 

wash-out pattern defined by I(t1) < I(t2) and is determined by ^^^^^^^^^^^^ -^-^ - ^^'^ 

pattern is for example coded with the color red. If the intensity function = ^(^1) ^(^Q) 

which is a measure of the apparent initial rate of wash-in, has a high value for a slow 
wash-out pattern, than the red color will be bright reflecting fast wash-in (Fig. 13). If the 
intensity function has a low value for a slow wash-out pattern than the red color will be 
dark reflecting slow wash-in (Fig. 13). 

[0040] The second pattern is the moderate wash-out pattern (Fig. 14) defined 
[0041] by I(t1) «l(t2) is determined by -0.1 < ^^^^^^^^^ pattern is. for 

example, coded with the color green. As in the red case, if the intensity function has a 
high value for this moderate pattern, then the green color will be bright (Fig. 14). If, 
however, the intensity function has a low value, the green color will be dark (Fig. 14). 
The third pattern (Fig. 15) is the fast washout pattern defined by I(t1 ) > I(t2) and is 

determined by ^^^^^"^^^^^ < 0.1 . This pattern is coded for example blue. Again if the 
7(^0 

intensity function is high, namely, wash-in apparent initial rate fast, the color blue will 
be bright (Fig. 15). If, however, the intensity function is low, the color blue will be dark 
(Fig. 15). Thus, for each pattern of wash-out coded by a color hue, there can be 
defined a wash-in rate which relates to the intensity function which determines color 
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intensity. The separation between different wash-in rates depends on the range of 
color intensities chosen. 

[0042] Referring again to Fig. 2, in block 20 the apparatus, in the final output for the 
selected tO, t1 and t2 and for the selected system parameters, assigns to each grid 
point that defines values of the variables of the system a color hue and a color 
intensity. The colorized and intensity normalized grid points collectively are termed a 
"calibration map". If the composite calibration map for the selected tO, t1 , t2 and for the 
selected system parameters is satisfactory as will be explained in detail hereafter, the 
program ends and sends the output to storage in block 2 (Fig. 1 ) or directly to block 3 
(Fig. 1). If the calibration map is not satisfactory, the apparatus makes an appropriate 
adjustment of the three time points, for example, by incrementing tl and t2 and inputs 
into block 

[0043] 13. Alternatively, the apparatus can increment new values for parameters of the 
system and/or of the measurement, as preselected, and input into block 12. Also, both 
time points and these parameters can be changed. The steps in the program are 
repeated using the original inputs of blocks 10 and 1 1 and going from blocks 14 to 20. 
This iteration (block 21 in Fig. 2) can be repeated until a satisfactory calibration map is 
obtained for a set of preselected times and parameters. 

[0044] In the definition of the intensity function in block 16 (Fig. 2) and the definition of 
the pattern of wash-out/color hue function C in blocks 22 and 23 (Fig. 3) it is assumed 
that I(tl)>l(t0). In cases where l(tl) is negative but the absolute value of /I(t1)/ is higher 
than I(t0) the definitions hold for the absolute values. The choices for coloring function 
C in block 23 (Fig. 3) are not limited to the specific examples of C < -0.1, -0.1 <_ C < 
0,1 and C>_0.1 , but any other fraction (such as 0.05 or 0.2) can be used to define C 
depending on the system, 

i. Referring now to Fig. 4, there is shown a flow diagram 
(steps or blocks) for setting the parameters for data 
collection and creating a calibration map. This flow 
diagram is written for a monitored parameter that varies 
with time as a function of two variables of the system 
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assigned here with the letters K and v. For each grid 
point in a 2 dimensional grid of K and v, a pixel of 
dimension of 0.01 units of K and 0.01 units of v is defined 
in steps or bloclcs 38, 40, 53 and 55. The program starts 
in block 30 and gets inputs of the time points, tO, t1 and 
t2, system and measurement parameters and the range 
of K and the range of v between their min. and max. 
values in block 31. The program starts from pixel (K min, 
V min) in block 32 to calculate I(t0), I(t1), I(t2) in block 33 
using an approximate or exact equation correlated to the 
system being monitored, as is known, that estimates how 
the parameter monitored with time l(t) depends on K and 
V, and on other system parameters. The determined or 
calculated l(tl) and I(t0) are used to calculate for each 
pixel Intensity (K, v) as shown in block 34, which 
represent wash-in Initial rate. The program is then 
searching whether the pixel has max. Intensity (blocks 35 
to 40) and proceeds through all the pixels in a loop mode 
returning to block 33 and going again through the steps 
34 to block 40 until it reaches the pixel with maximum K 
and maximum v. Through this loop, the pixel with max. 
intensity is identified and intensity is calculated for all 
pixels (K, v). Now, the program proceeds to calculate for 
each pixel starting from pixel (K min, v min) block 43 a 
normalized intensity, block 44 nonnalized relative to the 
max. intensity. The pixel with max. intensity is assigned 
with a maximum value for intensity N. N can be 1 , 2, 3 or 
any number such as, 8, 64, 256 (computer numbers), etc. 
depending on the demands of the system. Then, the 
program calculates the wash-out pattern for each pixel 
starting from pixel (K min, v min) until it reaches pixel (K 



13 



max, V max) and codes with color hue each pattern as 
shown in blocks 45 to 54. Now, all pixels have a color 
hue and a normalized color intensity. This produces in 
the output a calibration map of K, v in block 56 for the 
selected tO, t1 , t2 and system and measurement 
parameters. If the calibration map is not satisfactory e.g. 
excessively slanted toward one color hue, new time 
points, or new system or measurement parameter 
values, or all are adjusted in the direction to correct the 
calibration map and bring it to a more satisfactory 
balance from a color distribution standpoint. The program 
goes through all the steps in the flow diagram again 
using the new inputs until a satisfactory calibration map is 
obtained, which sets the selected time points and system 
parameters. What will -be a satisfactory calibration map 
depends on the system and will be apparent to one 
skilled in the art. For most systems a satisfactory map will 
have about a third of the pixels red, a third green and a 
third blue. 

[0045] A specific example of a flow diagram for setting the parameters for tracer 
modulated NM, termed also contrast enhanced MRI, is shown in Fig. 5. The MRI signal 
is the monitored parameter that is changing with time as a result of administrating a 
tracer, termed also contrast agent. The input parameter in step 6 includes the three 
time points tO, t1 12. These time points are obtained initially by experience in step 87 
and are subjected to an iterative process until the best three time points are obtained. 
Other inputs are the tracer-dose and the MRI parameters that define how the NM 
signal is recorded. Both the tracer dose and MRI parameters can be constant, or can 
be optimized by the iterative process in step 87. The input also includes 
pharmacokinetics parameters that define the tracer change with time in the blood, and 
max. and min. values for the two variables K and v that define pathophysiological 
characteristics of the system, namely, a subject body. The variable K defines 
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microvascular permeability which estimates the capacity of blood vessels to leak out 
the tracer. The variable v defines the fraction of extracellular volume which estimates 
the amount of free space in a tissue. The steps 62 to 85 in this flow diagram follow the 
steps 32-54 in Fig. 4, In this flow diagram the maximum intensity is assigned in step 74 
to have the value of 256. The output in step 86 consists of a calibration map of the two 
variables K and v ranging between K min, v min to K max, v max for a specific set of 
time points and the other inputs. Each pixel in this map with specific K, v values has a 
color hue and a color intensity. A satisfactory calibration map is defined by reaching a 
certain distribution of the colors or of the colors and color intensities. For example, a 
satisfactory map can be a map that divides the K-v plan or plane, or volume between 
the three colors to approximately three equal areas, namely, approximately a third of 
the pixels in the calibration map are red, a third are green and a third are blue. 
[0046] Shown in Fig. 16 is a typical calibration map according to the present invention. 
The map was created based on the equations of contrast enhancement as known in 
the art, for the variables microvascular pemieability K and fraction of extracellular 
volume V. The map is constructed as an expected pattern (e.g.. color and color 
intensity) for any three or more selected time points. More explicitly, and with reference 
to Fig. 16, a three-dimensional representation in a two-dimensional drawing is used. 
Fig. 16 shows the time points for humans using a gradient echo and a tracer dose of 
0.08 mmol/kg, with the time points: tO = 0, t1 = 4 and t2 = 12 min. These time points 
were selected in order to discriminate between e.g., fibroadenoma and carcinoma. The 
isotherms represent regions of the same intensity in each pattern, e.g., same initial rate 
of wash-in. One dimension is microvascular permeability K ranging between values of 
interest (for example, 0.00 min"'* to 0.3 min'\ The second dimension is fraction of 
extracellular volume v ranging between 0.1 to 1 and the third dimension, normalized 
intensity is actually the intensity of each color (for example, any value between 0 to 
256 intensities). This calibration map serves to determine optimal preset time points tO, 
t1 and t2 and other parameters such as dose. For different systems it is possible to 
select different optimal time points. The calibration map also serves to interpret the 
output of the processing means of the apparatus, which for the specific example of 
MRI is a 3TP image, defined subsequently in the description. It is clear that if the three 
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time points are chosen to be very close together the calibration map will show only a 
slow wash-out behavior, namely, the red pattern according to the above example will 
predominate. On the other hand, if the last point is chosen very far in time, the 
calibration map will be dominated by fast wash-out. namely, dominated by blue. The 
suitable 3 time points for a specific system are selected by having all three colors 
distributed in the most 

[0047] revealing way, namely, in the calibration map about one third of the area is 
occupied by each color. 

[0048] Flow Diagrams similar to the flow diagram in Fig. 5 exist for other specific 
systems. At the final output, a satisfactory calibration map of the variables for an 
optimal set of tO, t1 , t2 and other inputs is always obtained. 

[0049] The time points and parameters are set at block 1 (Fig. 1) and are used by the 
control and effect means 3 to control and effect the data collection (Fig. 1). Referring 
now to Fig. 6, there is shown a detailed block diagram for controlling and effecting data 
collection. The input in blocks 90, 91 can be directly obtained from the storage 2 (Fig. 
1) or from block 1 (Fig. 1). The measurements in blocks 92, 95 and 97, are performed 
by means specific for each system. Such means are known in each art. There must be 
control of times of tracer administration and of measure so that tO, t1 and t2 of the input 
are accurately controlled. The injection or administration of tracer in block 94 into the 
system can be performed in any known way. The timing of the tracer administration is 
fixed to start at tO and should end usually, but not necessarily before t1 , preferably 
close to to. The data collected in blocks 92, 95 and 97 are transferred to processing 
means 5 (Fig. 1 ) either directly from block 3 (Fig. 1 ) or from storage 4 (Fig. 1 ). 
[0050] A specific example of a block diagram of a control and effect means for 
controlling and effecting the data collection, as part of a modified MRI. apparatus or 
machine for tracer modulated MRI, is shown in Fig. 7. The control and effect means 
shown in the block diagram of Fig. 7 receives as an input in block 120 the time points 
to, t1 and t2 set by the selecting means in block 1 in Fig. 1 , and described in detail for 
this specific example in the flow diagram 1 of Fig. 4. The other inputs relate to the MRI 
parameters and to the dose of the tracer and the phamriacokinetic parameters of the 
tracer that is injected into the blood of the subject. The MRI parameters and the dose 
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of tracer are set by the selecting means in block 1 of Fig. 1 , and are described in detail 
for tracer modulated MRI In the flow diagram of Fig. 4. Next, the apparatus, in block 

121 in Fig. 7 records an 

[0051] image of a defined area or a defined volume in the system, namely a body, by 
means known in this art. The recording parameters are those set by the selecting 
means 1 in Fig. 1 and inputted into block 12 Fig. 7. The image is then stored in block 

122 and a tracer, termed also a contrast agent, is administered at time tO into the body 
in block 123 by any known way. This administration is timed to start at tO for a duration 
that ends preferably, but not necessarily, close to time point tO and before recording at 
time t1 . After the administration of tracer, the apparatus, in block 124. at time tl records 
an image of the same area or volume as was recorded in block 121 in the same body 
using the same NM parameters as in block 121. This image is stored in block 125. 
Next, at time t2, the apparatus, records another image of the same area or volume in 
the same body, using the same NIRI parameters as in block 121 , and then stores this 
image in the storage of block 127. 

[0052] Referring now to Fig. 8 there is shown a more detailed block diagram of 
processing means shown in block 5 of Fig. 1 . The inputs shown in block 100 of Fig. 8 
are the time points tO, tl and t2 set in block 1 of Fig. 1 and used in block 3 of Fig. 1 . 
Another input shown in block 101 of Fig. 8 includes the data collected in block 3 of Fig. 
1 . These data can be directly transferred from block 3 of Fig. 1 or from the storage in 
block 4 in Fig. 1 . The data in block 101 of Fig. 8 for each spatial unit are presented by 
three data values obtained at three different times. The first data value for each spatial 
unit, is termed I(t0) and is obtained before administration of the tracer. The tracer 
and/or third is administered at time point tO. Thus, the first data value is measured 
before time point tO but as close to this time point as possible. The second data value 
for each spatial unit is 

[0053] obtained at time point tl and is termed I(t1). The third data value for each spatial 
unit is obtained at time point t2 and is termed I(t2). Another input in block 101 in Fig. 8 
is max Intensity. The max Intensity value is part of the output of the selected means as 
shown in block 20 in Fig. 2, 



17 



[0054] From the data I(t0), I(t1) and the time points tO, t1 the normalized intensity is 
then determined in block 102 for each spatial unit. The normalized intensity is given by 

the equation [I{tl)'-I(tO)]xN 

[I(tO)](tl - tO)x{Max Intensity) 

[0055] I(t0), I(t1) and Max Intensity have been inputted in block 101. N is an integer 
number that can be 1,2,3 or 8, 64, 256 (computer numbers). N in block 102 in Fig. 8 is 
equal to N in block 44 in Fig. 4, in the same monitoring of a system. Next, the 
apparatus in block 104 colorizes each spatial unit according to the program shown in 
Fig. 3. The final output in block 106 will have each spatial unit assigned with a color 
hue and a color intensity. The color hue represents a wash-out pattern and the color 
intensity represents an initial rate of wash-in. The colored output in block 106 is fed to 
storage 6 in Fig. 1 or directly to the analysis means in block 7 of Fig. 1 . For each color 
and color intensity in each spatial unit in the output in block 106 of Fig. 8 there is the 
same color and color intensity in at least one grid point or location in the calibration 
map with defined values of variables of the system. Thus, the color/color intensity in 
each spatial unit obtained by the processing means 5 can be related to defined values 
of variables detemriined in the calibration map. 

[0056] The apparatus shown in Fig. 1 functions with normalized intensities. The 
normalization is performed by selecting means in block 1 in Fig. 1 in the process of 
obtaining a calibration map defined in detail previously. The selecting means defines 
max intensity and uses the value of this max intensity to normalize all intensities 
measured to this max intensity. Alternatively, it is possible, but not preferred, to 
separate the normalization in reference to max intensity, performed by the selecting 
means, from that performed by the processing means, by choosing the maximum 
intensity independently by the selecting means and by the processing means. 
However, this weakens the correlation between the calibration map and the final output 
of data processing shown in block 106 of Fig. 8 and in block 126 of Fig. 9 and in block 
166 of Fig. 10. 

[0057] Referring now to Fig. 9, there is shown a flow diagram of steps or blocks of the 
apparatus for processing collected data. This diagram deals with data recorded in two 
dimensions so that each data point is a pixel. The first step 1 10 in Fig. 9 is the start. 
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This is followed by input in step 111. The input consists of the selected time points tO, 
t1 and t2 set by the selected means in block 1 in Fig. 1 and used by the control and 
effect means in block 3 in Fig. 1. The spatial units n and m define, a pixel in a grid. For 
the first pixel n =1 and m =1 and then n goes from 1 to n and m goes from 1 to m in 
steps of 1. Another input is the collected data I(t0), I(t1) and I(t2) of each pixel. This 
input can be directly obtained from the control and effect means in block 3 in Fig. 1 or 
from the storage in block 4 of Fig. 1 . The input also includes the max intensity which is 
obtained in the output of the selecting means shown in block 20 in Fig. 2 or in step 56 
in the flow diagram of Fig. 4. Next, the apparatus, in steps 112 and 1 13 in Fig. 9 
selects the first pixel n=m=1 and in step 1 14 detemiines the normalized intensity in this 

pixel defined by the intensity function 

N 

[0058] and by the normalization in reference to the max intensity of 



Max Intensity 

[0059] where N= integer > 1 as defined above. The normalized intensity is a measure 
of the initial rate for the wash-in behavior. N in step 1 14 is equal to N in step 44 in Fig. 
4 in the same monitoring of a system. Next, this pixel (n = m = 1 ) is colorized in steps 
1 1 5 to 1 20 using the color function C according to the block diagram in Fig. 3 and 
according to steps 45 to 50 in Fig. 4. This determines wash-out pattern/ color hue of 
this pixel. If none of the possible patterns (3 patterns) occur, for example, as 1(t1) = 0 
this pixel is colored in black in step 121. In steps 122, 123, 124, 125 the next pixel is 
selected and a loop to step 114 that follows until step 120 determines for this next pixel 
the normalized intensity and the pattern/color hue. This loop is repeated for all pixels 
[0060] (n X m). In the output of step 126 in Fig. 9 each pixel has a wash-out pattern 
coded by a color and a wash-in initial rate coded by color intensity. The output is 
followed by step t27 which ends the flow diagram of Fig. 9. 

[0061] A specific example of a flow diagram of apparatus for processing data collected 
for tracer modulated MRI is shown in Fig. 10. The flow diagram starts with step 150. 
The next, step 151 , is the input of the selected time points tO, t1 , t2 set by the selecting 
means of the apparatus as shown for tracer modulated MRI in the apparatus depicted 
in the flow diagram of Figure 5 and is part of the output step 86 in Figure 5. These 
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three time points are also used by the control and effect means for data collection in 
the same monitoring of the body in step 151 of Fig, 10. IVIaximum intensity is also 
obtained from the output step 86 in Fig. 4 and is obtained in the same monitoring of the 
system, namely, the body. The input also includes the data collected by the control and 
effect means as shown for tracer modulated MRI in Fig. 7. These data are recorded 
images. There are three recorded images or three sets of recorded images of the 
same area or volume. The first image or set of images is recorded prior to tracer 
administration. The second image or set of images is recorded at time tl and the third 
image or set of images is recorded at time t2. Each pixel or voxel in the image has an 
MRI signal intensity which is changing with time after tracer administration. In the first 
image recorded prior to tracer administration, but close to the administration time, the 
intensity is termed limagei(n.m)- The pixel for which n=1 and m=1 is called the first pixel 
with intensity ljmagei(i,i)- There are n x m pixels where n goes from 1 to n and m goes 
from 1 to m. The intensity in each pixel (n.m) in the second image recorded at time 
point tl is termed Iimage2{n.m)- Similarly the Intensity in each pixel (n,m) of the third 
image recorded at time point t2 is termed Iimage3(n,m)- In steps 152 and 153 in Fig. 10 
the pixels are assigned starting from n=1 and m=1. In step 154 in Fig. 10 the 
normalized color intensity is determined for pixel (n,m) 

[0062] from the values of pixel (n,m) in limagei(n.m)and Iimage2{n.m) according to: 

[0063] normalized intensity= ^("■")"^"^g^KH.n,) ^ 256 

[^«n.g.i(«.«)](^i -^o) maxint ensity 

[0064] The normalized intensities are determined for all pixels by a loop in steps 162, 
[0065] 163, 164 and 165 until all pixels have been processed. 

(1) Next, in steps 155 to 160 the wash-out pattern/color is determined for 

each pixel (n,m) from the values Iimage2(n.m) and Iimage3(n,m) using the color 

function C which 
[0066] is defined for this specific example as: 

1. C = Color Function = ^^'^^g^^C"'-) " ^w(n.m) 

^ imagelin^m) 
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[0067] and a color is chosen as shown in block 23 in Fig. 3, and as shown in steps 155 
to 159 in Fig. 10. Pixels that are left uncolored through 155-160 as is the case when 
[00681 Iimage2(n,m) = 0 are asslgned black in step 161. The colorizing steps are looped by 
steps 162, 163, 164, 165 to include all pixels. The output in step 166 shows an n x m 
image in which each pixel (n,m) has a defined wash-out pattern/color hue and a 
defined wash-in rate/color intensity. This colored image is termed the three time points 
image or, in short, 3TP image. The color hue and color intensity is correlated to the 
calibration map and interpreted in terms of the values of microvascular permeability K 
and fraction of extracellular volume v, the two variables of the subject body in the 
calibration map of the output in step 86 in Fig. 5. 

[0069] Referring now to Fig. 1 1 , there is shown a block diagram of the structure of the 
analysis means 7 of Fig. 1. The input in block 200 consists of the calibration map 
obtained by the selecting means 1 and included in the output in block 20 in Fig. 2 or in 
the output in block 56 in Fig. 4. The other input in block 201 of Fig. 1 1 is the output of 
the processing means shown in block 126 in Fig. 9 or in block 106 in Fig. 8. The input 
in block 200 and the input in block 201 are from the same monitoring of the subject 
system. The analysis in block 202 consists of analysis of distribution of colors and of 
color intensities using for example histograms. The analysis consists also of a 
correlation between the calibration map and the color/color intensity of each spatial unit 
of the input in block 201 . This correlation estimates values of the variables of the 
calibration map for each spatial unit in the input of block 201 . Finally, the apparatus in 
block 203 outputs the distribution analysis and the correlation with the calibration map. 
[0070] A specific example of that portion of the apparatus including the means for 
analysis of the 3TP image obtained in a tracer modulated MRI is shown in the block 
diagram of Fig. 12. The input in block 210 is the 3TP image. The analysis in block 212 
of the 3TP image consists of analyzing color distribution and color intensity distribution, 
such as determining how many pixels are colored red with a certain intensity and 
making a similar determination for the other colors. Also, a separate analysis of 
intensity distribution and of color distribution can be performed. The part of the 
apparatus in block 21 1 provides means for inputting the calibration map obtained by 
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the selecting means as shown for tracer modulated NM in Fig. 5. This calibration map 
is for the same tO, t1 , t2 and other shared parameters as the 3TP image. Analysis of 
the 3TP image in terms of the two variables of the calibration map K and v is shown in 
block 213 of Fig. 12. This correlates color hue/color intensity in each pixel of the 3TP 
image to the values of the pathophysiological parameters K-microvascular permeability 
and v-fraction of extracellular volume determined by the color hue/color intensity in 
the-calibration map. Thus, the analysis is performed in terms of distributions of the two 
pathophysiological variables in the area or volume imaged. Finally, the apparatus in 
block 214 stores in a store or outputs as digital signals or displays in a display device 
like a monitor or is fed to a printer and a color print is obtained one or a plurality of 3TP 
images/and the corresponding correlated calibration maps. 
[0071] Another specific example relates to control and monitoring apparatus for an 
irrigation system. One of the most frequently used modern techniques to irrigate or 
water large areas in an efficient way is by drip irrigators. 

[0072] In the planning of such a system there are parameters such as the dimensions 
of the pipes, the extent and size of dripping holes and the pressure and timing of 
irrigation that can be adjusted according to the needs. These parameters will overall 
determine the rate of water dripping per unit area assigned here by the letter K. 
However, another parameter which will determine the efficiency of the irrigation is the 
water apparent diffusion constant in the ground, assigned here by the letter v. This 
diffusion rate or constant depends on the physical and chemical properties of the soil in 
the ground that the water passes through. Namely, in regions with light soil, such as 
sand, the apparent diffusion constant will be high while in regions with heavy soil it will 
be slow. Thus, the diffusion rate varies over the field needed to be irrigated. By the 
apparatus of the present invention it is possible to estimate K and v and then optimize 
the irrigation efficiency. 

[0073] As a preliminary matter one needs to measure the amount of water per unit 
weight of soil. There are several ways to determine water content. One for example is: 
weigh accurately an amount of soil just after digging it. .:Dry the soil completely and 
then measure again the weight The loss in weight is equivalent to the amount of water 
in this sample. 
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[0074J The samples can be taken with a spatial resolution that varies depending on the 
size of the field and on the accuracy needed to be reached. For example, for a very 
large field of tens or hundreds of square kms, it is reasonable to divide the field into 1 
km' units. Thus each pixel in the final image of irrigation constructed by the novel 
apparatus of the present invention and in particular by the processing means will 
reflect behavior per 1 km^. A sample of soil should be taken from the middle or any 
other defined location in this unit area of 1 km^. The size of the sample can vary but 
can be small of about 1 gram. The depth from which the sample should be taken can 
vary according to the needs. For example, if the growth of the plants to be irrigated 
depends on the amount of water at the level of the roots, then the sample should be 
taken from this level. It is also possible to use the apparatus in 3D and take samples 
from the same area but at varying depth. Samples should be taken from approximately 
the same place (the size of the sample is much smaller than the overall unit area). 
[0075] The apparatus of the invention performs as follows. The area to be tested has 
the dripping system ready for test and the positions of sampling are assigned. At a 
time, just before the start of the operation of the irrigation system, samples are taken 
from all the assigned positions. Then, at time point tO the operation of the irrigation 
system is initiated for a pre-set time which ends before time point t1 which is 
determined by the calibration map. After the irrigation is stopped, a second sample is 
taken for measurement from each assigned position at time point t1 . Finally, at the 
pre-set time point t2 a third sample from each assigned position is taken for 
measurement. The amount of water Is then determined in the samples. For each 
position the water content in the three samples taken at time point tO, and at time 
points t1 and t2 will change according to a wash-in rate and wash-out pattern and will 
be characterized by the color hue/color intensity code as developed according to the 
novel apparatus and method of the present invention. The wash-in and wash-out 
behaviors depend on 

[0076] K = rate of water dripping per unit area and on v = water apparent diffusion 
constant in the ground. For the same v the initial rate will increase with K, for the same 
K the initial rate will increase with v. The wash-out pattern will also depend on K and v. 
If V is faster than K the fast wash-out pattern (assigned blue) will predominate. If K and 
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V are of the same order, the green color, coding moderate wash-out, will predominate. 
If Q is lower than K a slow wash-out pattern, red will predominate. 
[0077] The distribution of K, v over all positions (each position is described by a pixel) 
is detemnined by correlating the coloring and color intensity of the pixels (each with the 
color hue and color intensity) in reference to the calibration map. 
[0078] The calibration map is providing wash-in intensity function and wash-out pattern 
in a K-v plane for tO, t1 , t2 as follows: a calculation for each pair of K, v of the amount 
of water accumulated at time points t1 and t2 is performed based on a model known to 
those skilled in the art. The range of K is chosen from 0 to the maximum level of the 
subject irrigation system while the range of v is from 0 to the diffusion constant of pure 
water or, when known, the highest water diffusion constant in the field to be irrigated. 
To obtain optimal resolution within the range of K and v, time points tO, t1 and t2 are 
chosen in such a way that the K-v plane will be divided between the three wash-out 
patterns/colors to approximately three equal areas. 

[0079] Once measurements have been made and the novel apparatus of the invention 
has performed its initial function, the apparatus can be adjusted to vary K (K can be 
varied since it depends on the irrigation system) in such a way that the irrigation in 
each defined area represented by a pixel will be the most efficient. For example, the 
moderate wash-out pattern (green) may be preferred in order to provide a constant 
amount of water over a defined time period. By increasing K it is possible to move from 
the blue region to the green for the same v. To assure the achievement of the final 
adjustment, it is possible to change K and then run the apparatus for the same three 
time points and positions in the subject field. 

[0080] Another use of the color/color intensity coded map of irrigation is the preparation 

of a plan of planting by adjusting the kind of plants or the density of planting to the 

quality of irrigation dictated by the kind of soil and the irrigation system. 

[0081] The new apparatus of the present invention can be utilized to test and modify 

the air condition planning, either for heating or for cooling or both. Described now in 

detail is novel apparatus for controlling and/or monitoring a heating process. 

[0082] For an air condition system (heating, cooling or both) built for a whole structure 

having within rooms or defined spaces, such as a house, a factory, an office building. 



24 



shopping mall or a complex of houses, it is important to design the system in such a 
way that each defined space will be conditioned efficiently and then, to verify the 
reality. Certain regions may be overheated while others can be over cold. The 
adjustment by the novel apparatus for such a structure can be done by modifying the 
amount of heat per unit area and unit time assigned, in this example, with the letter K. 
The other variable that will determine the heating capacity is the, rate of heat transfer 
per unit area to the surrounding environment due to imperfection in the isolation, in this 
example assigned the letter v. 

[0083] The assessment of the heat is performed by measuring the temperature by 
thermocouples or thermometers placed at any number of locations within each room or 
defined space. Each thermometer position will define a pixel position in a 2D or 3D 
plan of the subject system (structure). 

[0084] At a time point prior to tO. with the system turned off, readings of all temperature 
measuring devices in all positions are taken, determining temperature TO. Then at time 
point to the air condition system is turned on for a pre-set time, the end of which is 
before and close to time point t1 . At time point t1 , after the system has been turned off. 
second readings of all devices are taken determining temperatures T1 . Then, after a 
second preset time ending at time point t2 third readings are taken determining in each 
position temperature T2. The temperature changes between time points tO, t1 and t2, 
namely, the difference between temperatures TO, T1 and T2, for each position of 
measurement, will depend on the amount of conditioning per unit area per time (K) 
released in the room and on the amount of heat lost or gained from the outside through 
the walls by diffusion (v). The later parameter v can be negative or positive depending 
on the direction of flow of heat between the environment and the position where 
measurement occurs. These changes can be described by a wash-in pattern of air flow 
and wash-out pattern of air flow. The pattern of wash-out is described by a color hue. If 
the temperature T1 at time t1 is higher than the temperature T2 at time t2 the washout 
process is defined to be fast and is assigned blue. If T1 is smaller than T2 the wash out 
process is defined as slow and is assigned red. If the temperature will remain the same 
T1 = T2 (within a range predetermined by the apparatus) the wash-out process is 
defined to be moderate and is assigned the green color. The brightness of the colors 
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will depend on the initial rate defined by -^^ (usually but not necessarily tO = 

ro(n-ro) 

0). The Initial rate will also depend on K and on v. For the same v, the initial rate will 
increase with K. For the same K the initial rate will decrease with increasing v and will 
increase with decreasing v also to negative values. If the place is not well isolated and 
the surrounding is colder, then v will be high and the change in temperature will follow 
the pattern of fast wash-out. If the isolation is good (v is small and close to 0). the 
change in temperature will follow the pattern of moderate wash-out. If the surrounding 
is warmer and the room is not well isolated v will have a high negative value and the 
change in temperature will follow the slow wash-out pattern. Using the novel color 
hue/color intensity concept of the present invention it will be possible to identify places 
that are not well conditioned, e.g. heated, and are not well isolated and places that are 
over conditioned, e.g. overheated. 

[0085] If the air conditioning is heating, the temperature is determined by the amount of 
heat reaching the place where the temperature measuring device or element, e.g. 
thermometer, is placed and by the amount of heat that leaves this place as a result of 
heat loss to or heat gain from the surroundings. In certain cases defined by the size of 
the room and the distribution of the heat source we can assume that during the heating 
time the heat flow is relatively fast and equilibration in the room is rapidly achieved. 
Thus, the temperature will depend on the total amount of heat produced during the 
heating time period. K will therefore range between 0 and the maximum capacity of the 
heating system. The flow from or to the surroundings is determined by the same 
insulation, namely the same v but with opposite signs. Thus v will range between -v to 
+v with the actual value determined for example by the value with no insulation. 
[0086] The time points tO, t1 and t2 are chosen by using a calibration map constructed 
based on an equation known to those skilled in the art that relates the change in 
temperature with time to K as described above. Although, the particular example 
discussed related to heating, the novel apparatus can be used with the same program 
and means to evaluate cooling using absolute values for the changes in temperature 
that determine wash-out patterns and wash-in initial rates. 
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[0087] The apparatus of the present invention includes a computer system operating 
electronically, optically or both having a memory, a central processing unit, a display, 
an input device for generating device event signals and coacting therewith software for 
use with the computer. The software (in binary or related form) comprises a computer 
usable medium having computer readable program code thereon including the 
program logic for implementing the various flow charts and block diagrams described 
above. Since the details of computers are well known in the art and because persons 
skilled in the art have sufficient expertise and knowledge to be capable of 
implementing the flow charts and block diagrams, a detailed description of the specific 
hardware has been omitted as superfluous and unnecessary to a full and complete 
understanding and appreciation of the present invention as described above. Those 
skilled in the art will be able to make and use the apparatus and method of the present 
invention from the detailed description and teachings contained herein. 
[0088] Summarizing the 3TP method with respect to contrast enhance MRI, for a given 
breast lesion, changes in MRI signal intensity (SI) reflect changes in the concentration 
of the contrast agent. The concentration, in turn, is predominantly determined by two 
pathophysiological parameters that characterize malignant tumors and differentiate 
them from benign ones. These parameters are: the product: (blood vessel surface 
area) x (permeability) per unit volume; and the extracellular volume fraction (EVF) 
accessible to the contrast agent. From the breast MRI images, the 3TP algorithm 
detects the SI at each location, pixel-by-pixel, for one pre-contrast time point and two 
post-contrast time points (hence the name Three Time Point, or 3TP). The algorithm 
then codes the SI changes between the three time points using color intensity and 
color hue as follows: 

(1 ) Color Intensity codes the rate at which the SI changes between the first 
and second time points with a resolution of 256 intensities where dark 
colors signify slow change and bright colors signify rapid change. 

(2) Color Hue is a measure of contrast washout and is coded depending on 
the SI change between images recorded at the second and third time 
points. 
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[0089] The color coded 3TP images are related to pathophysiology via the 
mathematical model described above so that color hue and intensity are related to the 
product, (vessel surface area x permeability), and to the extravascular space (EVF), 
respectively, yielding a "calibration map", used for interpretation. Cancers typically 
show more bright red regions, reflecting the presence of increased vessel permeability 
and higher cell density. Benign tumors and normal breast conditions typically 
demonstrate greater areas of blue, indicative of the presence of lower cell density and 
thus higher extracellular volume with diminished vascular permeability. 

i) Figure 17 shows in view (a) that MLO mammographic projection of a 49 yo 
woman does not show 15 mm cancer which was palpable and diagnosed by 
in-office FNA as infiltrating ductal cancer; in view (b) that subtraction sagital 
MRI (0 minute (pre-contrast) image subtracted from 2 minute image shows 
ring-enhancing, spiculated malignancy in the anterior breast (arrow), and 
posteriorly, a small intramammary lymph note also enhances, but has a 
characteristic morphology, including fatty hilus (arrowhead); in view (c) that 
3TP parametric similar to that of 1 b, but with overlay of colored pixels 
superimposed on MRI image by 3TP software based on a physiological 
model (described in text) reveals bright red indication of high probability for 
malignancy and this lesion was prospectively given a score of 5 (highly 
suspicious for malignancy); and in view (d) shows a calibration map 
corresponding to the patient shown in la. The calibration map shows color 
hue (red, green, blue) based on differences in signal intensity between the 
second and third Images of the three image 3TP image set. Color intensity is 
based on the difference in SI between the first and second images of the 
image set. As described in the text, areas of high vessel permeability x 
surface area and low extravascular volume fraction (EVF), typical of 
malignancy, are coded as bright red. For optimal discrimination of benign and 
malignant lesions, the imaging parameters are chosen to approximately 
divide the calibration map into equal areas of red, green and blue. 
[0090] Accuracy of diagnosis relies on judicious selection of the three time points. This 
selection is made by adjusting time points so that, given the imaging parameters, the 
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red, green and blue pixels are equally distributed in this map. For the imaging 
parameters we employed, we found that the two post contrast imaging points that 
provided this optimal "calibration map" occurred at two and six minutes after contrast 
injection. Attempting interpretation using other time points will shift the distribution of 
red, green and blue pixels toward a higher sensitivity or specificity. By seeking those 
points in which red, green and blue pixels were equally distributed throughout the 
calibration map, it is hypothesized that the most accurate diagnosis could be obtained. 
The color convention used, is that cancer is now associated with red pixels, and benign 
conditions with blue pixels, contrary to any earlier convention. 
MR Imaging 

a) For a specific example of clinical testing, imaging was performed at 1 .5 Tesia 
(GE Medical Systems, Waukesha, Wl) using a phased array breast coil (MRI 
Devices, Waukesha, Wl). A three dimensional gradient echo acquisition was 
employed using parameters: TR= 15 msec; TE=4.2 msec; flip angle = 30 
degrees FOV = 16-18 cm; matrix = 256x256; NEX=1 .0, and slice thickness = 2.2- 
3.0 mm. Seven consecutive image sets of 56 slices (interpolated from 28 slices), 
were obtained over 14 minutes and 45 seconds. 

i) Gadodiamide (Omniscan - Nycomed Laboratories. Princeton, NJ) was 
injected three minutes after the beginning of the scan series, i.e. 1 minute 
after the start of the second scan sequence. Contrast was administered at 2 
cc/sec, followed by 15 cc of saline flush, also delivered at 2 cc/sec, using an 
automated pump (Medrad Corporation, Indiancis, PA). 

Image Interpretation and Data Analysis 

[0091] For the specific example of clinical testing, MR images were sent by File 
Transfer Protocol (ftp) from the laboratory performing clinical trials to the research 
laboratory for analysis. Interpretation of images involved visually examining each of 
the slices computed from imaging time points at 0, 2 and 6 minutes for a coherent 
group of pixels which could indicate a lesion. Prior experience with the 3TP method 
has shown that when a lesion has > 15% red pixels, it is likely malignant; if few red 
pixels (<10%) are present, the lesion is likely benign. Benign lesions typically show a 
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high fraction of blue pixels (>50%) and low color intensity. In this regard, note Figure 
18 which shows in view (a), a MLO mammographic projection showing an ovoid focal 
lesion along the inferior mid breast (arrow) in a 69 yo woman; pathological diagnosis, 
via excisional biopsy, was benign breast tissue; and in view (b) a 3TP sagital MRI 
parametric image showing virtually all dark blue pixels indicating a lesion with low 
values of vessel permeability x surface area and extravascular space - indicators of 
benignity; the prospective 3TP diagnosis was benign (suspicion level 2). 
[0092] For many lesions, diagnosis, simply based on these criteria, is conclusive. For 
visually indetemninate cases, a parametric 3TP image can be computed using the 
additional data points, 4 and 8 minutes, available as part of a seven image set. This 
technique shifts the calibration map towards one in which sensitivity is increased, at 
the price of specificity, with the idea that misdiagnosing a malignancy as benign is 
worse than misdiagnosing a benign lesion as malignant. With the three new time 
points, if the green pixels become red, the lesion can be diagnosed as malignant. 
[0093] Figure 19 shows in view (a) a MLO mammographic projection showing an area 
of ill-defined density and architectural distortion (arrow) in a 64 yo woman; diagnosis, 
via excisional biopsy, was infiltrating ductal cancer; in view (b) subtraction MRI (image 
at 0 minutes (pre-Gd contrast) subtracted from image at time 6 minutes) showing an 
irregular area of enhancement corresponding to the mammographic lesion; and in view 
(c) top three images being magnified views from adjacent slices of the central portion 
of the SIR parametric image, calculated using 0, 2 and 6 minute MR images. The 
results show a visually indeterminate number of red pixels. The bottom three images 
are these same image locations, but 3TP images recalculated using 0, 4 and 8 minute 
MR images. There is a shift toward an increasing number of red pixels, indicating 
malignancy to be more probable than benignancy. The prospective 3TP diagnosis was 
malignant, at suspicion level 4. 

[0094] The final diagnosis was graded using a scale slightly modified from the BIRAD 
scale used for mammographic interpretation in the United States: 

1 . 1 = very likely to be benign 

2. 2 = probably benign 

3. 3 = indetenninate (6 month f/u) 
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4. 4 = possibly malignant (biopsy) 

5. 5 = very likely to be malignant (biopsy) 

[0095] In the clinical test noted above, the 3TP method correctly diagnosed 27 of 31 
malignant (grade 4 or 5) and 31 of 37 benign lesions (grade 1 . 2 or 3). The ROC curve, 
based on the BIRAD-like classifications scheme is shown in Figure 20. The ROC cun/e 
of Figure 20 was derived from 3TP data in which lesion location was supplied to the 
researcher, but no other clinical information was supplied. The researcher was then 
asked to supply a number from 1 to 5 indicating probability for malignancy or benignity. 
[0096] The area under the ROC curve, Az, was 0.91 1 and the standard deviation in Az 
was 0.036. Only one lesion was graded as indetemiinate (grade 3) - pathology 
showed a benign intraductal papilloma. The results, classified by mammographic 
lesion, were further categorized as follows: 

i) Lesion Type Sensitivity Specificity 

11) All 68 lesions 87% 84% 

iii) 45 solid masses 96% 82% 

iv) Microcalcifications 63% 81% 

[0097] There were four false negative results: three lesions showing microcalcifications 
without mass and one small solid lesion. The three foci of microcalcifications were all 
pathologically diagnosed as ductal cancer in situ (DCIS) - two intermediate grade (8 
and 14 mm); and one low grade (four ducts). One example of a false negative result is 
shown in Figure 21 which shows in view (a) an optical close-up of mammographic 
magnification CC view - 66 yo woman with cluster of microcalcifications in the superior 
mid breast; pathological diagnosis was DCIS, low to intemiediate grade. In view (b) of 
Figure 21 is shown an optical close-up of subtraction sagital MRI image - 6 minutes 
after injection of Gadodiamide; small irregular focus of enhancement (arrow) 
corresponds to cluster of microcalcifications. Figure 21 view (c) shows an optical 
close-up of sagital 3TP image showing an area of predominantly bright green pixels 
thought to be benign (score = 2). A total of three out of eight DCIS lesions were 
misdiagnosed as benign by the 3TP method. An area of future investigation is to 
determine whether there is a specific pattern for DCIS that will allow increased 
accuracy of diagnosis. 
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[0098] Note that five other foci of DCIS were correctly diagnosed as malignant. The 
misdiagnosed solid lesion measured 5 mm in diameter and had a pathological 
diagnosis of invasive ductal cancer. In the same breast, two larger lesions, measuring 
8 and 14 mm were correctly diagnosed as malignant. 

[0099] Six false positive results were obtained: one 1 1 mm focus of fibrocystic change; 
one 3 mm focus of schlerosing adenosis; one 9 mm intraductal papilloma; one 10 mm 
focus of mixed pathology (fibrocystic change and fibroadenoma); one 3 mm intraductal 
papilloma and one 9 mm fibroadenoma. Of fourteen fibroadenomas, 12 were correctly 
diagnosed as benign, while two were thought to be a malignancies. Figure 22 shows in 
view (a) ultrasound examination of a mammographically occult palpable lesion in a 45 
yo woman showing a gently tabulated mass, without acoustic shadowing, typically of 
fibroadenoma, which was confirmed by ultrasound-guided FNA. In Figure 22 view (b) 
subtraction MIP MRI, 6 minutes after Gadodiamide injection, shows that the 
fibroadenoma confinmed by ultrasound is the largest of multiple enhancing smaller 
lesions. At the workstation, many of these smaller lesions showed an enhancement 
profile similar to that of the larger, palpable and sonographically confirmed lesion. In 
clinical practice evaluation of the multiple other enhancing lesions by manual 
placement of an ROI is impractical. Though internal septations are said to be an 
important MRI sign of fibroadenoma, they were not noted in this patient. In Figure 22 
view (c) 3TP parametric image showed predominantly central bright green and 
peripheral blue pixels, consistent with a benign lesion (score = 2). The 3TP parametric 
image showed no other suspicious lesions. 

[0100] There were 7 patients (12.5%) who inadvertently benefited from having breast 
MRI. In four patients, a second or even third focus of malignancy was detected that 
was unexpected by mammography and changed the surgical approach. An example is 
shown in Figure 23 in which view (a) of an optical close-up of mammography 
(MLO projection) in a 44 yo woman showed a spiculated mass (straight arrow). An 
unexpected, mammographically occult, 2"^ lesion was detected at the location 
indicated by the curved arrow, as a result of this patient volunteering for the 3TP 
clinical trial. Pathological diagnosis in both cases was infiltrating ductal cancer. As is 
shown in Figure 23 view (b) a 3TP parametric sagital plane image of lesion suspected 
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to be malignant by mammography. High predominance of bright red pixels indicates 
high value of the product of vascular permeability x surface area, and low 
extravascular volume fraction, indicating high probability of malignancy (score = 5). In 
Figure 23 view (c) a 3TP parametric image of a second adjacent sagital slice showed a 
second site very suspicious for malignancy (score 5). Both sites were biopsied 
confirming unsuspected multifocal malignancy. 

[0101] In one patient, where the radiologist, based on mammographic results, 
suspected a 1 cm tumor, MRI revealed that it was actually 4 cm in size, this was 
subsequently diagnosed as an infiltrating ductal cancer. 

[0102] In one patient, in whom a well marginated lesion with rapid contrast washout 
was seen, the surgeon declined imaging guidance in removal of the palpable lesion. 
However, when the pathological report described only benign breast tissue, the 
surgeon, based on the MRI data, was urged to repeat the ultrasound, which confirmed 
lack of excision. Image-guided excision was then performed, now showing a benign 
papilloma. 

[0103] Finally, in one patient, after a failed mammographically-guided localization of 
mammographically vague and sonographically occult lesion, MRI-guided needle 
localization was used to excise the lesion and arrive at the diagnosis - invasive ductal 
cancer. 

[0104] For 68 pathologically proven lesions, the 3TP method, in a heterogeneous 
population, achieved an overall sensitivity of 87% and a specificity of 84% for detection 
of malignancy. Of note is that when results from the 45 solid masses were reviewed, 
the 3TP method achieved a sensitivity and specificity of 96% and 82%, respectively. 
Only one 5 mm malignancy was misdiagnosed as benign in a patient in whom two 
other larger malignant (8 and 14 mm) and one benign lesion (29 mm post lumpectomy 
seroma) were correctly diagnosed. 

[0105] Although the invention has been described in detail, nevertheless changes and 
modifications which do not depart from the teachings of the present invention will be 
evident to those skilled in art. Such changes and modification are deemed to come 
within the purview of the present invention and the appended claims. 
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